Euler-Bessel and Euler-Fourier Transforms 



Robert Ghrist 

Departments of Mathematics and Electrical /Systems Engineering, 
University of Pennsylvania, Philadelphia PA, USA 

E-mail: ghrist@math.upenn.edu 



Michael Robinson 

Department of Mathematics, University of Pennsylvania, Philadelphia 
PA, USA 

E-mail: robim@math.upenn.edu 



Euler Transforms 



2 



Abstract. We consider a topological integral transform of Bessel 
(concentric isospectral sets) type and Fourier (hyperplane isospectral sets) 
type, using the Euler characteristic as a measure. These transforms convert 
constructible Z-valued functions to continuous R-valued functions over a 
vector space. Core contributions include: the definition of the topological 
Bessel transform; a relationship in terms of the logarithmic blowup of 
the topological Fourier transform; and a novel Morse index formula 
for the transforms. We then apply the theory to problems of target 
reconstruction from enumerative sensor data, including localization and 
shape discrimination. This last application utilizes an extension of 
spatially variant apodization (SVA) to mitigate sidelobe phenomena. 



AMS classification scheme numbers: 65R10,58C35 
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Integral transforms are inherently geometric and, as such, are invaluable 
to applications in reconstruction. The vast literature on integral geometry 
possesses hints that many such integral transforms are at heart topological: 
note the constant refrain of Euler characteristic throughout integral- 
geometric results such as Crof ton's Theorem or Hadwiger's Theorems (see, 
L9J). The parallel appearance of integral transforms in microlocal 
analysis — especially from the sheaf-theoretic literature [8] — confirms the 
role of topology in integral transforms. 

This paper considers particular integral transforms designed to extract 
geometric features from topological data. The key technical tool involved 
is EULER CALCULUS — a simple integration theory using the (geometric or 
o-minimal) Euler characteristic as a measure (or valuation, to be precise). 
We define two types of Euler characteristic integral transforms: one, a 
generalization of the Fourier transform; the other, a generalization of the 
Hankel or Bessel transform. These generalizations are denoted Euler- 
FOURIER and Euler-Bessel transforms, respectively. These transforms are 
novel, apart from the foreshadowing in the recent work on Euler integration 
10. 

The contributions of this paper are as follow: 

(i) Definitions of the Euler-Bessel and Euler-Fourier transforms on a 
normed (resp. inner-product) vector space; 

(ii) Index theorems for both transforms which concentrate the transforms 
onto sets of critical points, and thereby reveal Morse-theoretic 
connections; 

(iii) Applications of the Euler-Bessel transform to target localization and 
shape-discrimination problems; and 

(iv) An extension of spatially variant apodization (SVA) to Euler-Bessel 
transforms, with applications to sidelobe mitigation and to shape 
discrimination. 

2. Background: Euler calculus 

The Euler calculus is an integral calculus based on the topological Euler 
characteristic. The reader may find a simple explanation of Euler calculus 
in llj and more detailed treatments in ffT0l[14l|. For simplicity, we work in 
a fixed class of suitably " tame " or DEFINABLE sets and mappings. In brief, 
an O-MINIMAL STRUCTURE is a collection of so-called definable subsets of 
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Euclidean space satisfying certain axioms; a definable mapping between 
definable sets is a map whose graph is also a definable subset of the product: 
see |[13|. For concreteness, the reader may substitute "semialgebraic" or 
"piecewise-linear" or " globally subanalytic" for definable. All definable 
sets are finitely decomposable into open simplices in a manner that makes 
the Euler characteristic invariant. Given a finite partition of a 
definable set A into definable sets a a definably homeomorphic to open 
simplices, 

X (A):=^(-l) dim ^. (1) 

a 

This Euler characteristic also possesses a description in terms of alternating 
sums of (local) homology groups, yielding a topological invariance (up to 
homeomorphism for general definable spaces; up to homotopy for compact 
definable spaces). 

The Euler characteristic is additive: x(^4u£?) = x(A)+x(B) — x(APiB). Thus, 
one defines a scale-invariant " measure " d\ and an integral via characteristic 
functions: J \a&X := x(A) for A definable. The collection of functions 
from a definable space X to Z generated by finite (!) linear combinations 
of characteristic functions over compact definable sets is the group of 
CONSTRUCTIBLE functions, CF(X). The Euler INTEGRAL is the linear 
operator J x • dx ' CF(X) — ^ Z taking the characteristic function l a of an 
open fc-simplex a to (— l) k . By additivity of x, the integral is well-defined 
[[TOl H4|. There are several formulae for the computation (and numerical 
approximation) of integrals with respect to dx Q. 

The integral with respect to dx is well-defined for more general 
constructible functions taking values in a discrete subset of R; however, 
continuously- varying integrands are problematic. A recent extension of the 
Euler integral to R-valued definable functions uses a limiting process Q. 
Let Def (X) denote the definable functions from X to R (those whose graphs 
in X x R are definable sets). There is a pair of dual extensions of the Euler 
integral, [dx\ an d \dx~], defined as follows: 

/ h\dx\ = lim - / \nh\dx, [ h\dx\ = lim - / \nh]dx- (2) 
Jx n ^°° n Jx Jx n ^°° n Jx 

These limits exist and are well-defined, though not equal in general. The 
TRIANGULATION THEOREM for Def (X) lO states that to any h G Def (X), 
there is a definable triangulation (a definable bijection to a disjoint union 
of open affine simplices in some Euclidean space) on which h is affine on 
each open simplex. From this, one may reduce all questions about the 
integrals over Def (X) to questions of affine integrands over simplices, using 
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the additivity of the integral. Using this reduction technique, one proves the 
following computational formulae Q: 

Theorem 2.1. For h e Def {X), 

p poo 

/ h[d X \ = / x{h>s}- X {h<s}ds (3) 

JX Js=0 
p poo 

/ h\dx] = / x{h>s}- X {h<-s}ds. (4) 

JX Js=0 

This integral is coordinate-free, in the sense of being invariant under right- 
actions of homeomorphisms of X; however, the integral operators are not 
linear, nor, since — J h \_dx\ = f h \dx\ , are they even homogeneous with 
respect to negative coefficients. The compelling feature of the measures 
\_dx\ and \dx\ is their relation to stratified Morse theory Q. 

Let C C X denote the set of critical points of h. For arbitrary p e C, the 
CO-INDEX of p, T*(p), is defined as 

r(p)= lim x(B e (p)n{h<h(p)+e'}), (5) 

where B e (p) denotes the closed ball in X of radius e about p. The dual INDEX 
at p is given by 

Z,(p) = lim x (B e (p) n{h> h(p) - e 7 }) , (6) 
Theorem 2.2 (Theorem 4 of Q). Ifh is continuous and definable on X, then: 

[ h[_d X \= [ hTd X) [ h\dx\= [ hl*d X . (7) 
Jx Jx JX JX 

This has the effect of concentrating the measure [^xj on the critical points of 
the distribution. In the case of h a Morse function on an n-manifold M, for 
each critical point p £ C Xf X\p) = (-l) dimM -Mrt andX*(p) = (-1)^), where 
ji(p) is the Morse index of p. Thus, if h is a Morse function, then: 

/ h ldx\ = E (-l) n -" {p) h(p) [ h \dx\ = £ (-ir ip) Kp) (8) 
Jm peC(h) Jm P eC(h) 



3. Definition: Euler-Bessel and Euler-Fourier transforms 

There are a number of interesting integral transforms based on dx, including 
convolution and Radon-type transforms ||4lfTT||. We introduce two Euler 
integral transforms on vector spaces for use in signal processing problems. 
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For the Eulerian generalization of a Bessel transform, let V denote a finite- 
dimensional vector space with (definable, continuous) norm ||-||, and let 
B r {x) denote the compact ball of points {y : \\y — x\\ < r}. Recall that 
CF denotes compactly-supported definable integer-valued functions. 

Definition 3.1. For h e CF(V) define the BESSEL TRANSFORM of h via 



This transform Euler-integrates h over the concentric spheres at x of radius 
r, and Lebesgue-integrates these spherical Euler integrals with respect to 
r. For the Euclidean norm, these isospectral sets are round spheres. Given 
our convention that CF(V) consists of compactly supported functions, B : 
CF(V r ) — » Z is well-defined using standard o-minimal techniques (the Conic 
Theorem |I3|). 

3.2. Fourier 

There is a similar integral transform that is best thought of as a topological 
version of the Fourier transform. This is a global version of the microlocal 
Fourier-Sato transform on the sheaf CF(V r ) [8J. For this transform, an 
inner product on V must be specified. The Fourier transform takes as its 
argument a covector (Gf. 

Definition 3.2. For h e CF(V) define the FOURIER TRANSFORM of h in the 
direction (gV* via 



JO Jt'Hr) 

Example 3.3. For A a compact convex subset of R n and ||£|| = 1, (J 7 1a)(0 
equals the projected length of A along the £-axis. 

The Bessel transform can be seen as a Fourier transform of the log-blowup. 
This perspective leads to results like the following. 

Proposition 3.4. The Bessel transform along an asymptotic ray is the Fourier 
transform along the ray's direction: for h e CF(V) and x ^ e V, 




(9) 




(10) 




(11) 



where x* is the dual covector. 
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Proof. The isospectral sets restricted to the (compact) support of h converge 
in the limit and the scalings are identical. □ 

While the Fourier transform obviously measures a "width" associated to a 
constructible function, the geometric interpretation of the Bessel transform 
is more involved. The next section explores this geometric content via index 
theory. 



4. Computation: Index-theoretic formulae 

The principal results of this paper is are index formulae for the Euler-Bessel 
and Euler-Fourier transforms that reduce the integrals to critical values. 

Lemma 4.1. For A cV the closure of a open subset ofV, star-convex with respect 
to x e A, 

Bl A (x) = [ d x \_d X l (12) 

JdA 

where d x is the distance-to-x function d x : V —> R + . 

Proof Consider the logarithmic blowup taking V — {x} = x ]R + , with 
the second coordinate being \\y — x\\. The level sets of the R + coordinate of 
the blowup are precisely the isospectral sets of B(x). The induced height 
function d x : d A —> M + is well-defined on the unit tangent sphere of x, 
8a = aS 72-1 , since A is star-convex with respect to x. By definition, 

poo 

Bh{x) = / X (AndB r (x))ds. 
Jo 

For A star-convex and top-dimensional, A n dB r (x) is homeomorphic to 
dA n {d x > r}. By Equation ©, 

Bh(x)= X (dAn{d x >r})dr= [ d x [d X \. 

JO JdA 

□ 

This theorem is a manifestation of Stokes 7 Theorem: the integral of the 
distance over dA equals the integral of the 'derivative 7 of distance over A. 
For non-star-convex domains, it is necessary to break up the boundary into 
positively and negatively oriented pieces. These orientations implicate [^xj 
and \dx\ respectively. 



Euler Transforms 



8 



Theorem 4.2. For A c V the closure of a definable bounded open subset of V 
and x e V, decompose dA into OA = d+Au d~A, where d^A are the (closure 
of) subsets of dA on which the outward-pointing half spaces contain (for d~) or, 
respectively, do not contain (for d+) x. Then, 

Bl A (x) = [ d x ld x \ - I d x \dx] (13) 



d x Td x - / d x I*d x . (14) 

C x nd+A Jc x nd~A 

where C x denotes the critical points of d x : dA —> [0, oo). 

Proof Assume, for simplicity, that A is the closure of the difference of C+, 
the cone at x over A+, and C~ , the cone over A~. (The case of multiple 
cones follow by induction.) These cones, being star-convex with respect 
to x, admit analysis as per Lemma 14.11 The crucial observation is that, by 
additivity of x, 

X (dB r (x) nA)= X {dCt n {d x > r}) - x(dC~ n {d x > r}). 

Integrating both sides with respect to dr and invoking Theorem 12.11 gives 

Bl A {x) = d x [d\\ - I d x \dx\. 

By Theorem 12.21 this reduces to an integral over the critical sets of d x . 
The only critical point of d x on C+ — dA or C~ — dA is x itself, on which 
the integrand d x takes the value and does not contribute to the integral. 
Therefore the integrals over the cone boundaries may be restricted to d + A 
and d~ A respectively. The index-theoretic result follows from Theorem 
IZ21 □ 

In even dimensions, the ldx\-vs-\dx\ dichotomy dissolves: 

Corollary 4.3. For dim V even and A the closure of a bounded definable open set, 

Bl A {x) = d x [dx\ = / d x Z*dx- (15) 

JdA JC X 

Proof For dim V even, dim dA is odd. Equation [18] of [2 J implies that on an 
odd-dimensional manifold, f \dx\ = — f L^xJ • Equation ([13)) completes the 
proof. □ 



Given the index theorem for the Euler-Bessel transform, that for the Euler- 
Fourier is a trivial modification that generalizes Example 13.31 




Figure 1. The index formula for B [left] and T [right] applied to 
1a localizes the transform to (topological or smooth) tangencies of the 
isospectral sets with dA. 



Theorem 4.4. For A cV the closure of a definable bounded open subset ofV and 
(GT- {0}, decompose dA into OA = d+A U d~ A, where d^A are the (closure 
of) subsets ofdA on which £ points out of(d + ) or into (d~) A. Then, 

^ia(o = / e L^xj - / e \dx\ (i6) 

Jd+A ^ d i A 

= f trd x - f i%d X . (17) 

Jc^nd+A Jc^nd'A 

where Q denotes the critical points of ^ : dA —> [0, oo). For dimV even, this 
becomes: 

^1a(0= / Cldx\ = [ £Z>d X . (18) 

J8A JC{: 

The proof follows that of Theorem 14.21 and is an exercise. Figure [1] gives a 
simple example of the Bessel and Fourier index theorems in M 2 . 

By linearity of B and T over CF(V), one derives index formulae for 
integrands in CF(V) expressible as a linear combination of 1^. for A4 the 
closure of definable bounded open sets. For a set A which is not of 
dimension dim V , it is still possible to apply the index formula by means 
of a limiting process on compact tubular neighborhoods of A. 

The remainder of this paper explores applications of the Euler-Bessel 
transform to signal processing problems involving target detection, 
localization, and discrimination. 



5. Application: Target localization 

Applications of Euler-type integral transforms are naturally made in the 
context of (reasonably dense) sensor networks. As in [P, we consider 
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the setting in which a finite number of targets reside in a field of sensors 
whose locations are parameterized by a vector space V. Each target has 
a corresponding SUPPORT — the subset of V on which a sensor senses 
the target, albeit without information of the target's range, bearing, or 
identity. The resulting sensor counting function h e CF(V) contains highly 
redundant but informative structure. The Euler integral can be used to 
extract from h the number of targets [GDI; the Euler-Bessel transform can be 
used to localize the targets in the case of convex target supports. Assume 
for example that all target supports are equal to a round ball: all sensors 
within some fixed distance of the target detect it. The Bessel transform of 
the sensor field reveals the exact location of the target. 

Proposition 5.1. For A = B R (p) a compact ball about p e K 2? \ the Bessel 
transform 31 a is a nondecreasing function of the distance to p, having unique 
zero at p. 

Proof. Convexity of balls and Corollary 14.31 implies that 



which equals diamA = 2R for x A and is monotone in distance-to-p within 



Note that Proposition 15.11 fails in odd dimensions; the Bessel transform of 
a ball in M 2n+1 is constant, and B obscures all information. However, for 
even dimensions, Proposition 15.11 provides a basis for target localization. 
For targets with convex supports (regions detected by counting sensors), 
the local minima of the Euler-Bessel transform can reveal target locations: 
see Fig. Efleft] for an example. Note that in this example, not all local 
minima are target centers: interference creates ghost minima. However, 
given h e CF(V r ), the integral intyhdx determines the number of targets. 
This provides a guide as to how many of the deepest local minima to 
interrogate. 

There are significant limitations to superposition by linearity for this 
application. When targets are nearby or overlapping, their individual 
transforms will have overlapping sidelobes, which results in uncertainty 
when the transform is being used for localization. A typical example of this 
difficulty (present even in the setting of round targets) is shown in Figure 



Bl A (x) = I d x ld X \ 



= max d x — min d, 

OA OA 



JdA 



A. 



□ 



Bright]. 
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Figure 2. The Euler-Bessel transform of a collection of convex targets [left] 
has local minima at the target centers. However, too much interference 
between targets obscures target centers [right]. 

6. Spatially variant apodization and sidelobes 

Figure El reveals the prevalence of SIDELOBES in the application of the Euler- 
Bessel transform — regions of " energy leakage " in the transform — much 
the same as occurs in Lebesgue-theoretic integral transforms. These are, 
in general, disruptive, especially in the context of several targets, where 
multiple side lobe interferences can create "ghost" images. Such problems 
are prevalent in traditional radar image processing, and their wealth of 
available perspectives in mitigating unwanted sidelobes provides a base of 
intuition for dealing with the present context. 

One obvious method for modifying the output of Euler-Bessel transform is 
to change the isospectral contours of integration , as regulated by the norm 
|| • || . For example, one can contrast the circular (or £ 2 ) Euler-Bessel transform 
with its square (or i^) variant. The resulting outputs of norm- varied 
transforms can have incisive characteristics in some circumstances: see §0 
However, it may the case that no single norm is optimal for a given input, 
especially if it consists of multiple targets with different characteristics. 

We propose the adaptation and refinement of one tool of widespread 
use in traditional radar processing. This usually goes under the name 
of SVA: SPATIALLY VARIANT APODIZATION (see [12J for the original 
implementation; an updated discussion is in [5]). Though there are many 
heuristic implementations of SVA, the core concept behind the method 
involves using a parameterized family of kernels and optimizing the 
transform pointwise with respect to this family. The family of kernels is 
designed so as to reduce as much as possible the magnitude of the sidelobe 
phenomena, while preserving as much as possible the primary lobe. 
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Figure 3. Before (left) and after (right) application of SVA to the £oo Euler- 
Bessel transform of a pair of disjoint square targets. 



For the setting of the Euler-Bessel transform, we propose the following. 
Consider a parameterized family A of norms ||-|| a , a e A. The SVA Euler- 
Bessel transform is the pointwise infimum of transforms over A. 



where B r ^ a (x) is the radius r ball about x in the a-norm. For example, A 
could describe a cyclic family of rotated £ OQ norms. As shown in Figure [3l 
SVA with this rotated family eliminates the sidelobes from the transform of 
a square rotated by an unknown amount, while preserving the response of 
another nearby target. 

As in the case of traditional radar processing, there is benefit to SVA 
implementation even when the transform contour is not similar to the 
expected target shape. For instance, contours result in strong sidelobes 
in the transform of a hexagon (Figure IHleft]). Application of rotated SVA 
(Figure Ht right]) does not eliminate the sidelobes, but does dramatically 
reduce their magnitude. 

7. Application: Waveforms and shape discrimination 

Motivated by matched-filter processing, we apply Euler-Bessel transforms 
to perform geometric discrimination from enumerative data. Since 
the Euler-Bessel transform of a target has minimal sidelobes when the 
isospectral contours agree with the target shape, it can be used to construct 
a SHAPE FILTER. To test this, a simulation was run on a function h = ^ ■ 1^. 
for Ai a collection of polygonal domains in M 2 . The results are shown in 
Figures [5] through [SJ In these, A\ (upper left) is a hexagon; A 2 (upper right) 




(19) 
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Figure 4. Before (left) and after (right) application of SVA to the ioo 
Euler-Bessel transform of a hexagonal target support. Use of SVA greatly 
mitigates sidelobes. 




Figure 5. Circular (£2) Euler-Bessel transforms of a collection of interfering 
target supports reveals the location of the round target. 
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Figure 6. Square (£oo) Euler-Bessel transforms of a collection of interfering 
target supports reveal the location of the square target. 



is a rotated square; A 3 (lower left) is a round disk; and A4 (lower right) is an 
axis-aligned square. The sizes of these support regions were varied, to the 
point of significant overlap and interference. 
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Figure 7. SVA (rotated £oo) Euler-Bessel transforms of a collection 
of interfering target supports reveal the location of the square targets, 
independent of orientation. 




-100 -50 50 100 -100 -50 50 100 -100 -50 50 



Figure 8. When target overlaps are too great, false minima appear in the 
Euler-Bessel transforms {£2, £00, and SVA ^oo) 

Three variants of the Euler-Bessel transform were applied to the functions: 
an £ 2 transform, an transform , and an SVA (rotated ^oo) transform. 
The deep local minima of each transform correspond to the likely centers 
for the selected target shapes, in accordance with Proposition [57T1 For 
instance, with the £ 2 transform, although there are minima at the center 
of the two squares, they are not as deep as the minimum at the center 
of the disk. Similarly, the square transform has its deepest minimum 
at the center of the axis-aligned square, and the SVA transform detects 
both squares. Even when the support regions overlap, the Euler-Bessel 
transforms still successfully discriminate target geometries, though too 
much overlap generates interference and obscures the target geometries 
(Figure [8]). 

Future work will explore the use of SVA in the context of Euler integral 
transforms. 
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(i) The index formulae (Theorems 14.21 and I4.4|) for the Bessel and 
Fourier transforms act as a localization of the integral transform and 
lead to closed-form expressions; this is essential to the computa- 
tions/ simulations in the present paper. 

(ii) We have left unaddressed several important issues, such as the 
existence and nature of inverse and discrete transforms. For inverses 
to a broad but distinct class of Euler integral transforms, see [ITTl ?]. 

(iii) Computational issues for discrete Euler-Bessel and Euler-Fourier 
transforms remain a significant challenge. Though individual Euler 
integrals can be efficiently approximated on planar networks H), the 
integral transforms of this paper are, at present, computation-intensive. 

(iv) There is no fundamental obstruction to defining T and B over Def (V), 
so that the transforms act on real- valued (definable) functions. In this 
case, one would have to distinguish between \_dx\ an d \d\\ versions 
of the transforms. It is not clear what such transforms measure - 
geometrically or topologically - or which index formulae might persist. 
The passage to definable inputs means that these transforms will no 
longer be linear (as neither \_dx\ nor \dx\ is)- 

(v) The use of SVA methods is one of many possible points of contact 
between the radar signal processing community and Euler integration. 
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